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We propose an explanation for the electronic nematic state observed recently in parent iron-based 
superconductors [T.-M. Chuang et at, Science 327, 181 (2010)]. We argue that the quasi-one- 
dimensional nanostructure identified in the quasiparticle interference (QPI) is a consequence of 
the interplay of the magnetic (tt, 0) spin-density wave (SDW) order with the underlying electronic 
structure. We show that the evolution of the QPI peaks largely reflects quasiparticle scattering 
between electronic bands involved in the SDW formation. Because of the ellipticity of the electron 
pocket and the fact that only one of the electron pockets is involved in the SDW, the resulting 
QPI has a pronounced one-dimensional structure. We further predict that the QPI crosses over 
to two-dimensionality on an energy scale, set by the SDW gap, which we estimate from neutron 
scattering data to be around 90 meV. 
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Introduction. One of the key challenges in condensed 
matter physics is to understand the nature of the many- 
body states which manifest themselves in experimentally 
observable anomalous properties. Such states may ex- 
hibit subtle, or even entirely novel, forms of static or 
fluctuating order. One of the most recent examples 
is the nematic electronic structure observed by means 
of spectroscopic imaging-scanning tunneling microscopy 
(SI-STM) experiments in Ca(Fei_ x Coa : )2As2, parent ma- 
terial of iron-based superconductors [1]. Given a certain 
similarity between the phase diagrams of iron-based and 
cuprate superconductors [2| - both contain an antifcrro- 
magnetic phase at small, and a superconducting phase at 
larger, dopings - and in view of the checkerboard elec- 
tronic pattern observed earlier in the cupratesQ, this 
experiment also refocuses attention on a possible role of 
quasi-one-dimensional physics. Indeed, the possibility of 
nematic order arising from orbital physics in such quasi- 
two-dimensional electronic systems has already been dis- 
cussed in Ref.[4j. 

Despite such similarities, there are considerable differ- 
ences in the normal state electronic structure of the re- 
spective parent compounds. The iron-based compound 
exhibits two circular hole pockets of unequal size, cen- 
tered around the T point (0,0), and two elliptic elec- 
tron pockets centered at the (0, ±7r) and (±7r, 0) points 
of the unfolded Brillouin zone (UBZ, based on the Fe- 
lattice) 0-0I- Electron and hole bands are significantly 
nested, i.e. e£ ~ — e^+Q, where Qi is either Qi = (0,7r) 
or Q2 = (tt, 0). As nesting enhances SDW instabilities, 
several researchers have argued for an itinerant descrip- 
tion of the magnetism in those compounds based at least 
partly on nesting Q. Also, in this picture the specific se- 
lection of a (0, tt) or (tt, 0) magnetic order, as well as the 
anisotropy of the spin wave spectra were attributed to 
the ellipticity of the electronic pockets [l(| 



Here, we analyze signatures of this SDW order in SI- 
STM measurements like those reported in Ref. [l]. We do 
so within the framework of QPI developed in the context 
of the cuprates [l2|, EH • We show that the evolution of 
the QPI peaks is largely due to quasiparticle scattering 
between the electronic bands involved in the SDW for- 
mation. Because of the ellipticity of the electron pockets, 
and due to the fact that only one of the electron pockets 
is involved in the SDW, the resulting QPI pattern has a 
pronounced one-dimensional structure. This is in good 
agreement with the abovementioned experiments Our 
theory predicts a crossover to two-dimensionality in the 
QPI to occur at a scale set by (twice) the SDW gap, 2Ai, 
which we estimate to be around 90mcV. The exact value, 
however, depends on the size of the magnetic moments 
which can vary from compound to compound. 

This paper is organized as follows. We set the stage by 
introducing model, notation and parameters to describe 
the underlying electronic structure. We next provide a 
qualitative account of the gross features of SI-STM mea- 
surements based on the electronic structure of the SDW 
phase. This we then back up with a detailed calculation 
of the relevant Green functions within the T-matrix for- 
malism. We close with the discussion of the interplay of 
inter- and intra-band impurity scattering. 

The model. We employ an effective mean-field four 
band model with two circular hole pockets at (0, 0) (a- 
bands) and two elliptic electron pockets at Qi and Q2 
(£-bands)[Io|: 
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FIG. 1: (color online) Constant energy intensity 
maps of the spectral density, J^o- ImTr Go CT (k, uj), 
(left panel) and absolute value of the QPI, 
X) kCT ImTr [G 0CT (k,cj)i CT (k, k + q, w)Go -(k + q, u)] for non- 
magnetic (middle panel) and magnetic (right panel) 
impurities obtained as described in the text. The arrows de- 
note the characteristic scattering wave vectors which appear 
in the SDW state. The color bars refer to the intensity in 
units of states/eV. 



fit, e^ 2 = e + tp ([1 - e] cos(p x ) + [1 + e] cos(p y + tt)) - 
fi%. e accounts for the ellipticity of the electron pock- 
ets. Following our previous analysis of the spin wave 
excitations, we use Fermi velocities and size of the Fermi 
pockets based on Refs. namely t a = O.&beV, tp — 
-0.68ey, m = IMeV, )i 2 = l.UeV, e = OMeV, and 
e = 0.5. For these values, the Fermi velocities are 0.5eVa 
for the ai-band, where a is the Fe — Fe lattice spacing, 
and v x = 0.27eVa and v y — OAQeVa along x- and y- 
directions for the /3i-band, and vice versa for fc. We use 
a x = a y = a = 1. 

We introduce the experimentally observed (tt, 0) SDW 
order parameter, within a standard mean-field approx- 
imation: Ai oc I]p( c i lP «5 c /3ip+Qi70 ; <57}- in this state, 
one of the a fermions couples with only one band of j3 
fermions, leaving the other hole and electron bands - 
and hence their electron and hole FSs - unaffected by 
the SDW. Without loss of generality we direct Aj along 
the z- quantization axis. 

Results. In the following we focus on the discussion of 
the QPI in the SDW state at energies below twice the 
SDW gap. Throughout the paper we set Ai w 45meV, 
from a previous analysis of the experimental data[icij|. 

The simplest way to understand QPI qualitatively is 
to consider the evolution of the spectral function in the 



SDW state in the fashion pioneered for the cuprates[12j. 
In FigQJa), (d), (g) we show constant energy scans of 
the spectral density for positive energies. Zero energy, 
FigOJa), corresponds to the SDW-state Fermi surface 
(FS). Its C2 symmetry, lower than the symmetry of the 
normal state Fermi surface, is immediately apparent. Ba- 
sically, the FS consists of one hole pocket around the 
T-point, elliptic electron pockets at (±tt, 0) both not 
involved in the SDW, and two small electron pockets 
around the T-point that arise due to the folding of one 
hole and the other elliptic pocket at Qi. 

It then seems natural to expect that the anisotropy 
of the QPI will arise due to bands involved in the SDW 
formation, that is, by the inter- and intra-pocket scatter- 
ing shown by wave vectors qi - q3- This scattering will 
necessarily reflect the GVsymmetry induced by SDW or- 
der. In particular, q2 and q3 refer to intra-pocket scat- 
tering, while qi represents inter-pocket scattering from 
the edges of boomerangs which have the largest density 
of states (DOS). With increasing energy the scattering 
between the small pockets starts to interfere with that 
arising from the large electron pocket, FigOJd)-(g). The 
anisotropy of the spectral density induced by (0, tt) SDW 
order persists to larger energies. However, around ~ 2Ai, 
the influence of the SDW gap disappears (not shown) and 
the four-fold symmetry of the electronic structure is ef- 
fectively restored. 

The actual QPI which is believed to be measured in 
SI-STS[14( arises from quasiparticle scattering by pertur- 
bations internal to the sample such as non-magnetic or 
magnetic impurities. In order to back up the above qual- 
itative picture, we therefore perform a standard analysis 
of such processes based on a T-matrix description [l3j . In 
particular, we introduce the impurity term in the Hamil- 
tonian 
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where V^ k , and J™ a , define the non-magnetic and the 
magnetic point-like interaction term between the elec- 
trons in bands i and i' , respectively. In the following we 
orient the magnetic impurity in the z-direction (we did 
not find a dramatic change in our results by using a gen- 
eral orientation of the impurity). By changing k + Qi 
to k in c^jk+Qicr and defining the new Nambu spinor as 

= ( C a 2 kt' C Qikf> C /3ikf C Lkf C Q 2 ki' C mkU C Sikl' C LkJ.) 

we can write the Hamiltonian as 
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FIG. 2: (color online) Plots of QPI along AFM, q x , (a) and 
FM, q y , (b) directions. The bottom curve is at -50 meV and 
the top curve is at +50 meV. Consecutive curves are separated 
by 10 meV. Red, green and black curves are guides to the eye. 



where <g> is the direct product of matrices and 
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Here, we assume that the intraband impurity scattering, 
Mo, is bigger than the interband scattering between the 
bands separated by a large Q, uq, and set uq — 0.2u . 
The Green function matrix is obtained via G kk ' (t) = 
-(T^(t)^,(0)), whence 

G kk <K) = G k K)[<S kk , + WK)G k ,K)], (4) 

where G k (w„) = [iuj n — $kj is the bare Green's 
function of the conduction electrons. Solving the 
Dyson equation for the T- matrix £ kk '(w n ) = ?7 kk ' + 
5Zk" ^k k "G k //(w n )t k » k /(o; n ), the LDOS is obtained via 
analytic continuation ioj n — > E + i0 + according to 
N c (E,r) = -±ImTr[G(r,r,u;OL„^. B+ « ) +- > Note that 
interference between the two partial waves give rise to a 
spatial modulation of the amplitude of the total wave 
which is then reflected in the local density of states 
(LDOS). We now assume the Born approximation holds 
and set uq = 0.1t a . 

We show the absolute value of the resulting QPI for 
positive energies in FigQJb),(e),(h) for a non-magnetic 
impurity (7=1, 7'=0)and in FigQJc),(f),(i) for a mag- 
netic impurity (7=0, 7'=1), respectively. Overall the 
QPI map resembles well the structure anticipated from 
the spectral density maps. In particular, we find the 
overall G2 symmetry of the QPI maps whose structures 
are determined by scattering at momenta shown on the 
left panel. In particular, at lOmeV the QPI shows peaks 
at q2 and q3 which refer to intra-pocket scattering, as 
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FIG. 3: (color online) Real space image of the QPI at -7meV 
away from the impurity site. The inset shows the correspond- 
ing QPI image. The QPI interference originating from the 
structure in the spectral density is reflected in the strong 
~ 10 -T- 15a modulation along the FM [x] -direction and a weak 
~ 5 -7- 7a modulation along the AF (j/)-direction. In order to 
model the real Co impurity we adopt the impurity potential 
as a mixture of the magnetic (20%) and non-magnetic( 
parts. Intensity refers to states/eV. 



well as qi which denote the inter-pocket scattering orig- 
inating from the edges of the boomerangs. Overall, the 
QPI shown in Figjljb), (e) and Fig |3] resembles the one 
found experimentally in Ref. [1]. In particular, we find 
that as a result of SDW induced electronic structure re- 
construction, the scattering interference modulations are 
strongly unidirectional and show G2 symmetry with qi 
and q3 induced structure along the AF q x axis. In the 
ferromagnetic q y axis modulations with wave vector q2 
are also found. As in the electronic structure, this per- 
sists only up to a scale set by the size the SDW gap. 

To analyze the energy dispersion of the QPI interfer- 
ence we show in FigJ5] its evolution along q x and q y - 
directions for energies between -50 and 50 meV. The 
peaks associated with q2 along the ferromagnetic q x di- 
rection and qi and q3 along the AF q y direction are dis- 
persive and their velocities are in direct correspondence 
with the quasiparticle group velocities of the bands in- 
volved in the SDW formation. 

Finally, in FigJ3] we show a real space image of the 
QPI away from the impurity position and at -7meV. We 
observe that the QPI modulation is reflected in periodic 
structures seen along the x and y directions. In the AF y 
direction with qi structure, we find a weak 'stripe' pat- 
tern with a periodicity of the order of ~ 7a. In addition, 
we also find a stronger modulation of 10 — 15a along the 
FM x direction, consistent with the structure described 
as nematic in Ref. [l[ . The x — y asymmetry of the struc- 
ture is again a consequence of the SDW with (0, tt) mag- 
netic order. We also stress that these extra modulations 
arise purely due to scattering between the bands involved 
in the SDW - no extra folding of the bands occurs. 

One natural question to ask is: How stable is the ob- 
served QPI and the resulting real space 'stripe' structure? 



4 




FIG. 4: (color online) Constant energy intensity maps of 
the spectral density, ^ a ImTr Go<r(k, ui), (left panel) and 
QPI, E^ImTr[G? 0o -(k,w)t -(k ) k + cbw)Go -(k + q,a;)] for 
non-magnetic (middle panel) and magnetic (right panel) im- 
purities and negative energies. The arrows denote the char- 
acteristic scattering wave vectors which appear in the SDW 
state. The color bars refer to the intensity in units of 
states/eV. 

First, the QPI is sensitive to the size of the SDW gap. 
For large enough Ai the FSs of the bands involved in the 
SDW are completely gapped and the same is also true for 
the spectral densities at low energies. In this case, QPI 
will be determined by the bands which are not involved 
in the SDW, and, therefore, its structure will not show 
strong quasi-lD character. This may explain why the 
pure AF SDW state does not show any ^-symmetric 
structure in the parent compounds where the magnetic 
moment (and the corresponding SDW gap) is quite large. 
Only when it is reduced upon doping, the bands involved 
in the SDW are located close to the FS and the QPI 
structure, described above, becomes visible. Note that 
the SDW gap used in our calculations corresponds to 
/i « 0A[1b- An additional factor why doping might be 
crucial for an observation of the superstructure is that 
the impurities act as pinning centers and the interfer- 
ence is more pronounced once disorder is increased. We 
also note that as soon as SDW forms due to the nesting of 
electron and hole bands, the main role of the underlying 
orbital structure is to modify the absolute intensities of 
the QPI as this depends on the underlying orbital matrix 
elements. However, it is unlikely that these factors would 
restore the four-fold symmetry in the QPI. Instead, it is 
expected that the corresponding QPI will be the same as 
in our Fig. 1 but with extra intensity modulation due to 
orbital matrix elements 15]. 

We note in passing that an increase of the intraband 
impurity scattering, uq, does not change the results sig- 
nificantly in any way except that at large enough val- 
ues of uq, one finds local resonances around the impurity 
site. In addition, the change of the ratio between the in- 
terband and intraband scattering is more subtle although 



the resulting QPI still shows the C2 anisotropy even up to 
uq = uq. However, there is still another interesting effect 
that we find upon changing the energy from positive to 
negative values. In particular, in Figf?]we show the spec- 
tral density and the corresponding QPI from — lOmeV to 
— 60meV. At low negative energies, i.e. above -20meV, 
the structure remains the same as for OmeV. However, 
for -30meV we find a rotation of the QPI by 90 degrees. 
This occurs due to particle-hole asymmetry of the bands 
contributing to the SDW. In the SDW state, new energies 

are = i (V + ± \/& 1 -4 1 + Qi) 2 +* A l) ■ 

While E c produces small pockets along x (AF) direc- 
tion which are pronounced for positive energies; E d yields 
pockets along the y (FM) direction visible at relatively 
large negative energies. As a consequence, the QPI ro- 
tates 90 degrees at energies lower than -30meV. This pre- 
diction would be interesting to check experimentally. 

In summary, we have presented a theory for SI-STM 
measurements in the iron-based superconductor parent 
compounds in the presence of (0, it) SDW magnetic order. 
We find that QPI introduced by scalar non-magnetic, as 
well as magnetic, impurities gives rise to a periodic mod- 
ulation of the real space LDOS. Because of the ellipticity 
of the electron pockets, and the fact that only one of the 
electron pockets is involved in the SDW, the resulting 
QPI has a pronounced quasi-one-dimensional structure in 
good agreement with recent SI-STM experiments [lj]. We 
further predict that the QPI becomes two-dimensional at 
energies larger than twice the SDW gap, 2Ai, which as 
we argue from the analysis of the neutron scattering data 
should be ~ 90mcV. 
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